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Abstract. We introduce a constraint-based framework for studying in- 
finite qualitative simulations concerned with contingencies such as time, 
space, shape, size, abstracted into a finite set of qualitative relations. To 
define the simulations we combine constraints that formalize the back- 
ground knowledge concerned with qualitative reasoning with appropriate 
inter-state constraints that are formulated using linear temporal logic. 
We implemented this approach in a constraint programming system 
(ECL B PS e ) by drawing on the ideas from bounded model checking. The 
implementation became realistic only after several rounds of optimiza- 
tions and experimentation with various heuristics. 

The resulting system allows us to test and modify the problem spec- 
ifications in a straightforward way and to combine various knowledge 
aspects. To demonstrate the expressiveness and simplicity of this ap- 
proach we discuss in detail two examples: a navigation problem and a 
simulation of juggling. 



1 Introduction 
1.1 Background 

Qualitative reasoning was introduced in AI to abstract from numeric quan- 
tities, such as the precise time of an event, or the location or trajectory of an 
object in space, and to reason instead on the level of appropriate abstractions. 
Two different forms of qualitative reasoning were studied. The first one is con- 
cerned with reasoning about continuous change in physical systems, monitoring 
streams of observations and simulating behaviours, to name a few applications. 
The main techniques used are qualitative differential equations, constraint prop- 
agation and discrete state graphs. For a thorough introduction see |15| . 

The second form of qualitative reasoning focuses on the study of contingen- 
cies such as time, space, shape, size, directions, through an abstraction of the 
quantitative information into a finite set of qualitative relations. One then relies 
on complete knowledge about the interrelationship between these qualitative re- 
lations. This approach is exemplified by temporal reasoning due to pp, spatial 
reasoning introduced in |10| and |19j . reasoning about cardinal directions (such 
as North, Northwest); see, e.g., J7j, etc. 



In this paper we study the second form of qualitative reasoning. Our aim is to 
show how infinite qualitative simulations can be naturally formalized by means 
of temporal logic and constraint satisfaction problems. Our approach allows us to 
use generic constraint programming systems rather than specialized qualitative 
reasoning systems. By a qualitative simulation we mean a reasoning about 
possible evolutions in time of models capturing qualitative information. One 
assumes that time is discrete and that only changes adhering to some desired 
format occur at each stage. Qualitative simulation in the first framework is 
discussed in |16) . while qualitative spatial simulation is considered in |SJ. 

1.2 Approach 

In the traditional constraint-based approach to qualitative reasoning the qualita- 
tive relations (for example overlap) are represented as constraints over variables 
with infinite domains (for example closed subsets of 1Z 2 ) and path-consistency 
is used as the constraint propagation; see, e.g., 

In our approach we represent qualitative relations as variables. This allows us 
to trade path-consistency for hyper-arc consistency which is directly available in 
most constraint programming systems, and to combine in a simple way various 
theories constituting the background knowledge. In turn, the domain specific 
knowledge about simulations is formulated using the linear temporal logic. These 
temporal formulas are subsequently translated into constraints. 

Standard techniques of constraint programming combined with techniques 
from bounded model checking can then be used to generate simulations. To 
support this claim, we implemented this approach in the constraint programming 
system ECL l PS e . However, this approach became realistic only after fine-tuning 
of the translation of temporal formulas to constraints and a judicious choice of 
branching strategy and constraint propagation. To show its usefulness we discuss 
in detail two case studies. In each of them the solutions were successfully found 
by our implementation, though for different problems different heuristic had to 
be used. 

The program is easy to use and to interact with. In fact, in some of the case 
studies we found by analyzing the generated solutions that the specifications were 
incomplete. In each case, thanks to the fact that the domain specific knowledge is 
formulated using temporal logic formulas, we could add the missing specifications 
in a straightforward way. 

1.3 Structure of the paper 

In Section[21we discuss examples of qualitative reasoning and in Section[3]explain 
our formalization of the qualitative reasoning by means of constraints. Next, in 
Section 0] we deal with qualitative simulations by introducing inter-state con- 
straints which connect different stages of simulation and determine which sce- 
narios are allowed. These constraints are defined using linear temporal logic. 
Their semantics is defined employing the concept of a cyclic path borrowed from 



2 



the bounded model checking approach (see [5]) for testing validity of temporal 
formulas. 

In Section[S]we explain how the inter-state constraints are translated to con- 
straints of the underlying background knowledge. Next, in Section [S] we discuss 
technical issues pertaining to our implementation that generates infinite quali- 
tative simulations. In the subsequent two sections we report on our case studies. 
Finally, in Section 03 we discuss the related work. 



2 Qualitative Reasoning: Setup and Examples 

As already said, in qualitative reasoning, one abstracts from the numeric quan- 
tities and reasons instead on the level of their abstractions. These abstractions 
are provided in the form of a finite set of qualitative relations, which should 
be contrasted with the infinite set of possibilities available at the numeric level. 
After determining the 'background knowledge' about these qualitative relations 
we can derive conclusions on an abstract level that would be difficult to achieve 
on the numeric level. The following three examples illustrate the matters. 

Example 1 (Region Connection Calculus). The qualitative spatial reasoning with 
topology introduced in [TiJ\ and JU| is concerned with the following set of qual- 
itative relations: 

RCC8 := {disjoint, meet, overlap, equal, covers, contains, covered-by, inside}. 

The objects under consideration are here spatial regions, and each region pair is 
in precisely one RCC8 relation; see Fig. ^ 

The background knowledge in 

this Case is the Set Of possible rela- covered-by inside 

tion triples pertaining to triples of _^ 
regions. For example, the relation 
triple (meet, meet, meet) is possible 
since there exist three regions pair- 
wise touching each other. In contrast, dis j° mt covers contains 
the triple (inside, inside, disjoint) is 

impossible since for any three regions Fig. 1. The eight RCC8 relations 
A, B, C, if A is inside B and B is 
inside C, then A cannot be disjoint 

with C. The set of possible triples is called the composition table; it is pre- 
sented in the above two papers. In total, the table lists 193 relation triples. 

Example 2 ( Cardinal Directions) . Qualitative reasoning dealing with relative di- 
rectional information about point objects can be formalized using the set of 
cardinal directions 

Dir := {N,NE,E,SE,S,SW,W,NW, EQ}, 
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that consists of the wind rose directions together with the identity relation de- 
noted by EQ; see ^21- The composition table for this form of qualitative reasoning 
is provided in |17) . 

Example 3 (Relative Size). Qualitative reasoning about relative size of objects 
is captured by the relations in the set 

Size := {<, =, >}. 

The corresponding composition table is given in 

Other examples of qualitative reasoning deal with shape, directional infor- 
mation about regions or cyclic ordering of orientations. In some of them the 
qualitative relations are non-binary and the background knowledge is more com- 
plex than the composition table. To simplify the exposition we assume in the 
following binary qualitative relations. 

3 Formalization of Qualitative Reasoning 

In what follows we follow the standard terminology of constraint programming. 
So by a constraint on a sequence x±, . . . , x m of variables with respective do- 
mains dom(xi), . . . , dom(x m ) we mean a subset of dom(xi) x • • • x dom{x m ). 
A constraint satisfaction problem (CSP) consists of a finite sequence of 
variables X with respective domains and a finite set of constraints, each on a 
subsequence of X. A solution to a CSP is an assignment of values to its variables 
from their domains that satisfies all constraints. 

We study here CSPs with finite domains and solve them using a top-down 
search interleaved with constraint propagation. In our implementation we use 
a heuristics-controlled domain partitioning as the branching strategy and 
hyper-arc consistency of JS| as the constraint propagation. 

We formalize the qualitative reasoning within the CSP framework as follows. 
We assume a finite set of objects 0, a finite set of binary qualitative relations Q 
and a ternary relation CT representing the composition table. Each qualitative 
relation between objects is modelled as a constraint variable the domain of which 
is a subset of Q. We stipulate such a relation variable for each ordered pair of 
objects and organize these variables in an array Rel which we call a qualitative 
array . 

For each triple a, b, c of elements of O we have then a ternary constraint comp 
on the corresponding variables: 

comp(Rel[a, b], Rel[b, c], Rel[a,c\) := 

CT n (dom(Rel[a,b]) x dom(Rel[b, c]) x dom(Rel[a,c])). 

To assume internal integrity of this approach we also adopt for each ordered 
pair a, 6 of elements of O, the binary constraint conv(Rel[a, b], Rel[b, a]) that 
represents the converse relation table, and postulate that Rel[a,a] = equal for 
all aeO. 

We call these constraints integrity constraints. 
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4 Specifying Simulations using Temporal Logic 




In our framework we assume a conceptual neighbourhood between the qual- 
itative relations. This is a binary relation neighbour between the elements 
of the relation set Q describing which atomic changes in the qualitative re- 
lations are admissible. So only 'smooth' transitions are allowed. For exam- 
ple, in the case of the Region Connection Calculus from Example the 
relation between two regions can change from disjoint to overlap only indi- 
rectly via meet. The neighbourhood relation for RCC8 has 22 elements such 
as (disjoint, meet), (meet, meet), (meet, overlap) and their converses and is shown 
in Fig. 

We assume here that objects can 
change size during the simulation. If wc 
wish to disallow this possibility, then the 
pairs (equal, covered-by), (equal, covers), 
(equal, inside), (equal, contains) and their 

converses should be excluded from the M t - L ^ 

conceptual neighbourhood relation. disjoint contains 

In what follows we represent each covers 
stage t of a simulation by a CSP T> t 

uniquely determined by a qualitative ar- pig _ 2 Thg RCC8 neighbourhood 
ray Qt and its integrity constraints. Here re i a tion 
t is a variable ranging over the set of natu- 
ral numbers that represents discrete time. 

Instead of Qt[a, b] we also write Q[a, b, t], as in fact we deal with a ternary array. 

The stages are linked by inter-state constraints that determine which sce- 
narios are allowed. The inter-state constraints always include constraints stipu- 
lating that the atomic changes respect the conceptual neighbourhood relation. 
Other inter-state constraints are problem dependent. 

A qualitative simulation corresponds then to a CSP consisting of stages all of 
which satisfy the integrity constraints and the problem dependent constraints, 
and such that the inter-state constraints are satisfied. To describe the inter-state 
constraints we use atomic formulas of the form 

Q[a, b] e U, Q[a, b] <£ TZ, Q[a, b] = q, Q[a, b] ± q, 

where TZ C Q and q € Q. As the latter three forms reduce to the first one, we 
deal with the first form only. 

We employ a propositional linear temporal logic with four temporal op- 
erators, O (eventually), O (next time), □ (from now on) and U (until), and 
with the usual connectives. We use bounded quantification as abbreviations, 
e.g., 4>{o\) V ... V 4>(ok) abbreviates to 3 A e {oi, ... ,0k}- 4>{A). 

Given a finite set of temporal formulas formalizing the inter-state constraints 
we wish then to exhibit a simulation in the form of an infinite sequence of 'atomic' 
transitions which satisfies these formulas and respects the integrity constraints. 
In the Section [3] we explain how each temporal formula is translated into a 
sequence of constraints. 
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Qi Q2 Qt Qk 

Fig. 3. A {k - £)-loop 

Paths and loops We now proceed by explaining the meaning of a temporal 
formula <f> with respect to an arbitrary infinite sequence of qualitative arrays, 

7r := Qx, Q2, ■ ■ ■ , 

that we call a path. Our goal is to implement this semantics, so we proceed in 
two stages: 

— First we provide a definition with respect to an arbitrary path. 

— Then we limit our attention to specific types of paths, which are unfoldings 
of a loop. 

In effect, we use here the approach employed in bounded model checking; see [5]- 
Additionally, to implement this approach in a simple way, we use a recursive 
definition of meaning of the temporal operators instead of the inductive one. 

We write \=n </> to express that <j> holds along the path it. We say then that 7r 
satisfies cj>. Given ir := Qx, Q2, ... we denote by 7Tj the subpath Qi, Qi+i, ■ ■ ■■ 
Hence wx = tt. The semantics is defined in the standard way, with the exception 
that the atomic formulas refer to qualitative arrays. The semantics of connec- 
tives is defined independently of the temporal aspect of the formula. For other 
formulas we proceed by recursion as follows: 



Ki Q[a,b] en 


if 


Q[a, b, i] G 1Z; 




if 


Ki+i 0; 


K 00 


if 


j=7T (j) and (=,r OHip; 


K 04> 


if 


\=„ or \=„ OO0; 


Kxu^ 


if 


K <t> or K xAO(xU 



Next, we limit our attention to paths that are loops. Following 5 we call a 
path 7r := Qx, Q2, ... a (A; — £)-loop if 

7r = u ■ v* with u := Qx, ■ ■ ■ , Qe-x and v :— Qe, . . . , Qk', 

see Fig. |3] By a general result, see [5], for every temporal formula if a path 
exists that satisfies it, then a loop path exists that satisfies <j>. This is exploited 
by our algorithm. Given a finite set of temporal formulas ^ it tries to find a path 
7r := Qx, Q2, ■ ■ ■ consisting of qualitative arrays that satisfies all formulas in 
by repeatedly trying to construct an infinite (k — £)-loop. Each such (k — £)-\oop 
can be finitely represented using k qualitative arrays. The algorithm is discussed 
in Section El 
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5 Temporal Formulas as Constraints 

A temporal formula restricts the sequence of qualitative arrays at consecutive 
stages (time instances). We now show how to translate these formulas to con- 
straints in a generic target constraint language. The translation is based on un- 
ravelling the temporal operators into primitive Boolean constraints and primitive 
constraints accessing the qualitative arrays. Furthermore, we discuss a variation 
of this translation that retains more structure of the formula, using non-Boolean 
array constraints. 

We assume that the target constraint language has primitive Boolean con- 
straints and reified versions of simple comparison and arithmetic constraints. 
(Recall that a reified constraint generalizes its base constraint by associating 
with it a Boolean variable reflecting its truth.) 

Paths with and without loops. Both finite and infinite paths can be accommo- 
dated within one constraint model. To this end, we view a finite sequence of 
qualitative arrays together with their integrity constraints as a single CSP. The 
sequence Qi, ■ ■ ■ ,Qk can represent both 

an infinite path tt = (<Qi, . . . , Qe-i) ■ (Qe, ■ ■ ■ , Qk)*, for some £ ^ 1 and k > £, 
or a finite path n = Qi, . . . , Qk- 

To distinguish between these cases, we interpret £ as a constraint variable. We 
define t = fc+1 to mean that there is no loop, so we have dom{£) — {1, . . . , fc+1}. 
A new placeholder array Qk+i is appended to the sequence of qualitative arrays, 
without integrity constraints except the neighbourhood constraints connecting 
it to Qk- Finally, possible looping is realized by conditional equality constraints 

(l = j) - (Qj = Qk+i) 

for all j E {1, . . . , fc}. Here Q p — Q q is an equality between qualitative arrays, 
i.e., the conjunction of equalities between the corresponding array elements. 

Translation into constraints. We denote by cons((j>,i) = b the sequence of con- 
straints representing the fact that formula 4> has the truth value b on the path 
7T,. The translation of a formula on Qi, . . . , Qk is initiated with cons(<j>, 1) = f . 
We define the constraint translation inductively as follows. 

Atomic formulas: 

cons(true,i) = b translates to b = 1; 

cons(Q[a\,a2\ £ R, i) = b translates to Q[a\, 02, i] = q,(q G R) = b. 
Connectives: 

cons(^(f>, i) = b translates to = b, cons(<p, i) = b'; 

other connectives are translated analogously. 
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Formula 0<f>: The next-time operator takes potential loops into account. 

cons(0(f>, i) = b translates to 

if i < k then 

cons((j), i + 1) = b; 
if i = k then 

£ = k + 1 -> b = 0, 

£ ^ k -> 6 = Aje{i,...,fe} = i -> cons(<t>,j)) . 

Formula O^: We translate 0</> by unravelling its recursive definition <fi\/00(fi. 
It suffices to do so a finite number n unravo i of steps beyond the current state, 
namely the number of steps to reach the loop, max(0, £ — i), plus the length 
of the loop, k — £. A subsequent unravelling step is unneeded as it would 
reach an already visited state. We find 

^unravel = k~ mm(£,i). 

This equation is a simplification in that £ is assumed constant. For a variable 
£, we 'pessimistically' replace £ here by the least value in its domain, min(£). 

Formulas D<f> and (f> U tp: These formulas are processed analogously to 0<fi. 

The result of translating a formula is a set of primitive reified Boolean con- 
straints and accesses to the qualitative arrays at certain times. 

Translation using array constraints. Unravelling the temporal operators leads 
to a creation of several identical copies of subformulas. In the case of the O tem- 
poral operator where the subformulas in essence are connected disjunctively, we 
can do better by translating differently. The idea is to push disjunctive informa- 
tion inside the variable domains. We use array constraints, which treat array 
lookups such as x = A[yi, . . . , y n ] as a constraint on the variables x, y\, . . . , y n 
and the (possibly variable) elements of the array A. Array constraints generalize 
the classic element constraint. 

Since we introduce new constraint variables when translating 0<fi using array 
constraints, one needs to be careful when 0<f> occurs in the scope of a negation. 
Constraint variables are implicitly existentially quantified, therefore negation 
cannot be implemented by a simple inversion of truth values. We address this 
difficulty by first transforming a formula into a negation normal form, using the 
standard equivalences of propositional and temporal logic. 

The constraint translations using array constraints (where different from 
above) follow. The crucial difference to the unravelling translation is that here i 
is a constraint variable. 

Formula 0</>: A fresh variable j ranging over state indices is introduced, mark- 
ing the state at which <f> is examined. The first possible state is the current 
position or the loop start, whichever is earlier. Both £ and i are constraint 
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variables, therefore their least possible values min(£), min(i), respectively, 
are considered. 

cons(0<j), i) = b translates to new j with dom(j) = {1, . . . , fc}, 

min(mwi(i), min(i)) ^ j, 
cons((p,j) = b. 

Formula 0</>: This case is equivalent to the previous translation of Ocf>, but 
we now need to treat i as a variable. So both "if . . . then" and — > are now 
implemented by Boolean constraints. 

6 Implementation 

Given a qualitative simulation problem formalized by means of integrity con- 
straints and inter-state constraints formulated as temporal formulas, our pro- 
gram generates a solution if one exists or reports a failure. During its execution 
a sequence of CSPs is repeatedly constructed, starting with a single CSP that is 
repeatedly step-wise extended. The number of steps that need to be considered 
to conclude failure depends on the temporal formulas and is finite [2]. The se- 
quence of CSPs can be viewed as a single finite CSP consisting of finite domain 
variables and constraints of a standard type and thus is each time solvable by 
generic constraint programming systems. The top-down search is implemented 
by means of a regular backtrack search algorithm based on a variable domain 
splitting and combined with constraint propagation. 

The variable domain splitting is controlled by domain-specific heuristics if 
available. We make use of the specialized reasoning techniques due to [201 f° r 
RCC8 and due to J7| for the cardinal directions. In these studies maximal 
tractable subclasses of the respective calculi are identified and corresponding- 
polynomial decision procedures for non-temporal qualitative problems are dis- 
cussed. In our terminology, if the domain of each relation variable in a qualita- 
tive array belongs to a certain class, then a certain sequence of domain split- 
tings intertwined with constraint propagation finds a solving instantiation of the 
variables without backtracking if one exists. However, here we deal with a more 
complex set-up: sequences of qualitative arrays together with arbitrary temporal 
constraints connecting them. These techniques can then still serve as heuristics. 
We use them in our implementation to split the variable domains in such a 
way that one of the subdomains belongs to a maximal tractable subclass of the 
respective calculus. 

We implemented the algorithm and both translations of temporal formulas to 
constraints in the ECI/PS e constraint programming system j^lj. The resulting 
program is about 2000 lines of code. We used as constraint propagation hyper- 
arc consistency algorithms directly available in ECI/PS e in its fd and propia 
libraries and for array constraints through the implementation discussed in [HJ. 
In the translations of the temporal formulas, following the insight from bounded 
model checking, redundancy in the resulting generation of constraints is reduced 
by sharing subformulas. 
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7 Case Study 1: Navigation 



Consider a ship and three buoys forming a triangle. The problem is to generate 
a cyclic route of the ship around the buoys. We reason qualitatively with the 
cardinal directions of Example |21 

— First, we postulate that all objects occupy different positions: 

□ Va,6GO. a^b -> Q[a,b] ^ EQ. 

— Without loss of generality we assume that the buoy positions are given by 
DQ[buoy al buoy c ] = NW, DQ[buoy a , buoy b ] = SW, DQ[buoy b , buoy c ] = NW 

and assume that the initial position of the ship is south of buoy c: 
Q[ship, buoy c ] = S. 

— To ensure that the ship follows the required path around the buoys we stip- 
ulate: 

□ (Q[ship 7 buoy c ] = S -» 0(Q[ship, buoy a ] = W A 

0(Q[ship, buoy b ] = N A 

O {Q[sbip, buoy c ] = E A 

0(Q[shi P , buoy c ]=S ))))). 

In this way we enforce an infinite circling of the ship around the buoys. 

When fed with the above constraints our 
program generated the infinite path formed ^-^ ^-^ ^-^ 

by the cycle through thirteen positions de- /v' ~* w "* vA, 

picted in Fig.^] The positions required to be _ / A 
visited are marked bv bold circles. Each of v!/ :" "\/ 



visited are marked by bold circles. Each of 
them can be reached from the previous one 
through an atomic change in one or more 



t 
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qualitative relations between the ship and ' 
the buoys. One hour running time was not ' © J'*" 
enough to succeed with the generic first-fail x s : 

heuristic, but it took only 20 s to find the © *~ © * 

cycle using the Dir-specific heuristic. The 

array constraint translation reduced this . XT . . , 

,. , , , _ Fig. 4. Navigation path 

slightly to 15 s. 

The cycle found is a shortest cycle satisfying the specifications. Note that 
other, longer cycles exist as well. For example, when starting in position 1 the 
ship can first move to an 'intermediate' position between positions 1 and 2, 
characterized by: 

Q[ship, buoy c ] = SW, Q[ship, buoy a ] = SE, Q[ship, buoy b ] = SE. 

We also examined a variant of this problem in which two ships are required 
to circle around the buoys while remaining in the N or NW relation w.r.t. each 
other. In this case the shortest cycle consisted of fifteen positions. 
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8 Case Study 2: Simulating of Juggling 

Next, we consider a qualitative formalization of juggling. We view it as a pro- 
cess having an initialization phase followed by a proper juggling phase which is 
repeated. As such it fits well our qualitative simulation framework. 

We consider two kinds of objects: the hands and the balls. For the sake of 
simplicity, we only distinguish the qualitative relations 'together', between a 
ball and a hand that holds it or between two touching balls, and 'apart'. This 
allows us to view the juggling domain as an instance of an existing topological 
framework: we identify 'together' and 'apart' with the relations meet and disjoint 
of the RCC8 calculus. 

In our concrete study, we assume a single juggler (with two hands) and three 
balls. We aim to recreate the three-ball-cascade, see ^] p. 8]. So we have five 
objects: 

O := Hands U Balls, 
Hands := {left-hand, right-hand} , 
Balls := { balk \ i E {1,2,3}}. 

The constraints are as follows. 

— We only represent the relations of being 'together' or 'apart': 

\J\/x,y E O. (x ^ y — ► Q[x,y] G {meet, disjoint}). 

— The hands are always apart: 

□ Q[left-hand, right-hand] = disjoint. 

— A ball is never in both hands at the same time: 

□ V6 S Balls. -> (Q[Mt-hand, b] = meet A Q [right-hand, b] = meet) . 

— From some state onwards, at any time instance at most one ball is in any 
hand: 

OD (V6 e Balls. V7i G Hands. Q[b, h] = meet 

V6 2 G Balls. b^b 2 — > V7i 2 € Hands. Q[6 2; h 2 ] = disjoint). 

— Two balls touch if and only if they are in the same hand: 

□ (V&i,& 2 £ Balls. h^b 2 -> 

(Q[foi,&2] = meet «-> 3h G Hands. (Q[h,bi] = meet A Q[h,b 2 ] = meet))). 

— A ball thrown from one hand remains in the air until it lands in the other 
hand: 

□ (V6 G Balls. Vfti, h 2 G Hands, hi hi A Q[h\, b] = meet — > 

Q[h\, b] = meet U (Q[hi, b] — disjoint A Q[h 2 , b] = disjoint A 
(Q[hi, b] = disjoint U Q[h 2 , b] = meet))). 
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— A ball in the air will land before any other ball that is currently in a hand, 

□ (V/ii, h 2 G Hands. V&i, b 2 G Balls. Q[hi, b\] = disjoint A Q[h 2 , b 2 ] = meet — » 
Q[h 2 , b 2 ] = meet U ((V/i G Hands. Q[h, b 2 ] = disjoint) 

U (3h G Hands. Q[h,bi] = meet))). 

— No two balls are thrown at the same time: 

□ (V&i, b 2 G Balls. b\ ^ b 2 — > Wii, /12 G Hands. 

->(Q[bi, hi] = meet A O Q[h, hi] = disjoint A 
Q[b 2 ,h 2 ] = meet A O Q[b 2 , h 2 ] = disjoint)). 

— A hand can interact with only one ball at a time: 

□ V/i G Hands. V&i G Balls. 

(Q[h, bi] = meet A O 61] = disjoint V 
Q[h, b%] = disjoint A O Q[/i, 61] = meet) — > 
V6 2 G Ba/Zs. 61 7^ 6 2 -> 

(Q[/i, 6 2 ] = meet -> O Q[/i, b 2 ] = meet) A 
(<9[/i, 62] = disjoint — » O 62] = disjoint). 

Initially balls 1 and 2 are in the left hand, while ball 3 is in the right hand: 

Q[left-hand, balh] = meet, Q[lcft-hand, ball 2 ] = meet, Q[right-hand, balls] = meet. 

Note that the constraints enforce that the juggling continues forever. Our 
program finds an infinite simulation in the form of a path [1..2] [3.. 8]*; see Fig- EI 
The running time was roughly 100 s using the generic first-fail heuristic; the 
RCC8-specific heuristic, resulting in 43 min, was not useful. 

We stress the fact that the complete specification of this problem is not 
straightforward. In fact, the interaction with our program revealed that the 
initial specification was incomplete. This led us to the introduction of the last 
constraint. 
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Fig. 5. Simulation of Juggling 
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Aspect Integration: Adding Cardinal Directions 

The compositional nature of the 'relations as variables' approach makes it easy 
to integrate several spatial aspects (e.g., topology and size, direction, shape etc.) 
in one model. For the non-temporal case, we argued in \J\ that the background 
knowledge on linking different aspects can be viewed as just another integrity 
constraint. Here we show that also qualitative simulation and aspect integration 
combine easily, by extending the juggling example with the cardinal directions. 

As the subject of this paper is modelling and solving, not the actual inter- 
aspect background knowledge, we only explain the integration of the three rela- 
tions meet, disjoint, equal with the cardinal directions Dir. We simply add 

link(Q[a, b], Qow[a,b\) := (Q[a, b] = equal) <-> (Qo\ r [a, b] = EQ) 

as the aspect linking constraint, ft refers to the two respective qualitative arrays 
and is stated for all spatial objects a, b. We add the following domain-specific 
requirements to our specification of juggling: 

Qo ir [lcft-hand, right-hand] = W; 

□ V& G Balls. Vft G Hands. Q[b, h] = meet -> Q Dir [b, h] = N; 

□ V6 £ Balls. Q[b, left-hand] = disjoint A Q [b, right-hand] = disjoint — > 

Q D ir[&, left-hand] ^ N A Q D ir[&, right-hand] ^ N. 

We state thus that a ball in a hand is 'above' that hand, and that a ball is not 
thrown straight upwards. 

This simple augmentation of the juggling domain with directions yields the 
same first simulation as in the single-aspect case, but now with the RCC8 and 
Dir components. The ball/hand relation just alternates between N and NW (or 
NE). 

We emphasize that it was straightforward to extend our implementation to 
achieve the integration of two aspects. The constraint propagation for the link 
constraints is achieved by the same generic hyper-arc consistency algorithm used 
for the single-aspect integrity constraints. This is in contrast to the 'relations 
as constraints' approach which requires new aspect integration algorithms; see, 
e.g., the bipath-consistency algorithm of |13| . 

9 Conclusions 

Related Work The most common approach to qualitative simulation is the 
one discussed in |151 Chapter 5] . For a recent overview see . It is based on a 
qualitative differential equation model (QDE) in which one abstracts from the 
usual differential equations by reasoning about a finite set of symbolic values 
(called landmark values). The resulting algorithm, called QSIM, constructs the 
tree of possible evolutions by repeatedly constructing the successor states. Dur- 
ing this process CSPs are generated and solved. This approach is best suited to 
simulate the evolution of physical systems. 
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Our approach is inspired by the qualitative spatial simulation studied in 9 , 
the main features of which are captured by the composition table and the neigh- 
bourhood relation discussed in Example Q The distinction between the integrity 
and inter-state constraints is introduced there; however, the latter only link con- 
secutive states in the simulation. As a result, our case studies are beyond their 
reach. Our experience with our program moreover suggests that the algorithm 
of [5] may not be a realistic basis for an efficient qualitative reasoning system. 

To our knowledge the '(qualitative) relations as variables' approach to mod- 
elling qualitative reasoning was first used in |21| . to deal with the qualitative 
temporal reasoning due to £Q. In |20| this approach is used in an argument to 
establish the quality of a generator of random scenarios, whilst the main part of 
this paper uses the customary 'relations as constraints' approach. In 2J pages 
30-33] we applied the 'relations as variables' approach to model a qualitative 
spatial reasoning problem. In [7] we used it to deal in a simple way with aspect 
integration and in [3] to study qualitative planning problems. 

In jS] various semantics for a programming language that combines temporal 
logic operators with constraint logic programming are studied. Finally, in the 
TLPlan system of 0] temporal logic is used to support the construction of 
control rules that guide plan search. The planning system is based on an incre- 
mental forward-search, so the temporal formulas are unfolded one step at a time, 
in contrast to the translation into constraints in our constraint-based system. 

Summary We introduced a constraint-based framework for describing infinite 
qualitative simulations. Simulations are formalized by means of inter-state con- 
straints that are defined using linear temporal logic. This results in a high de- 
gree of expressiveness. These constraints are translated into a generic target 
constraint language. The qualitative relations are represented as domains of con- 
straint variables. This makes the considered CSPs finite, allows one to use hyper- 
arc consistency as constraint propagation, and to integrate various knowledge 
aspects in a straightforward way by simply adding linking constraints. 

We implemented this approach in a generic constraint programming system, 
ECL l PS e , using techniques from bounded model checking and by experimenting 
with various heuristics. The resulting system is conceptually simple and easy to 
use and allows for a straightforward modification of the problem specifications. 
We substantiated these claims by means of two detailed case studies. 
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